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Images of flux vortices in superconductors acquired by transmission electron microscopy should 
allow a quantitative determination of their magnetic structure but so far, only visual comparisons 
have been made between experimental images and simulations. Here, we make a quantitative 
comparison between Fresnel images and simulations based on the modified London equation to 
investigate the magnetic structure of flux vortices in MgB2. This technique gives an absolute, low- 
field 30 Oe) measurement of the penetration depth from images of single vortices. We found that 
these simulations gave a good fit to the experimental images and that if all the other parameters 
in the fit were known, the penetration depth for individual vortices could be measured with an 
accuracy of ±5 nm. Averaging over 17 vortices gave a penetration depth of Kab ~ 113 ± 2 nm at 
10.8 K assuming that the entire thickness of the sample was superconducting. The main uncertainty 
in this measurement was the proportion of the specimen which was superconducting. Allowing for 
a non-superconducting layer of up to 50 nm thickness on the specimen surfaces gave a penetration 
depth in the range Aat = 100-115 nm, close to values of 90 ± 2 nm obtained by small-angle neutron 
scattering and 118-138 nm obtained by radio-frequency measurements. We also discuss the use of 
the transport of intensity equation which should, in principle, give a model-independent measure of 
the magnetic structure of flux vortices. 

PACS numbers: 74.25. Uv, 68.37.Lp 
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I. INTRODUCTION 



Superconductors expel magnetic flux from their inte- 
riors (the Mcissner effect) but if a magnetic field is ap- 
plied to a type-II superconductor which is larger than 
the lower critical field, Hd , magnetic flux penetrates the 
superconductor by flowing along channels called flux vor- 
tices. Each vortex carries a single quantum of magnetic 
flux given by $o = h/2e where h is Planck's constant and 
e the electron charge. They consist of a core where su- 
perconductivity is suppressed with a radius given by the 
coherence length, ^, surrounded by circulating supercur- 
rents which persist over a radius given by the penetration 
depth, A. When flux vortices move, energy is lost and so 
the performance of almost all superconducting devices is 
determined by the behaviour of flux vortices. 

Flux vortices can be imaged using transmission elec- 
tron microscopji^ as shown in Fig.[lja). The supercon- 
ductor is thinned to about 250 nm so that it is electron- 
transparent and the flux vortices penetrate normal to the 
thin surface. It is mounted at a high angle, a (typically 
45°), to provide a component of the B-fleld normal to 
the electron beam so that the electrons arc deflected by 
the Lorentz force and appear as black-white features in 
an out-of-focus image. Electron microscopy is the only 
imaging technique which gives information on the inter- 
nal magnetic structure of flux vortices, not just the stray 
flelds, and it allows imaging at video rate. As we show 
here, it can be used to investigate the magnetic structure 



of individual flux vortices and give an absolute measure 
of the penetration depth at low magnetic fields (typically 
30 Oe) irrespective of whether the vortices form a regular 
array or not. Obtaining an absolute value for the pene- 
tration depth is important as it gives information on the 
number density of electrons involved in superconductiv- 
itjs^, the nature of the superconducting stated and the 
types of vortex interaction which can occur—. 

The ability to measure the magnetic properties of indi- 
vidual vortices is useful for the investigation of materials 
in which the vortices do not form a regular array: an 
effect which is hampering current research into some of 
the iron-based superconductors^. The technique would 
also provide a good method to test predictions that vor- 
tices containing non-integer multiples of the flux quan- 
tum should occur in 2-componcnt superconductors such 
as MgB2 when their size becomes comparable to the co- 
herence lengthii^. It is also useful if the structure of 
certain vortices is altered by pinning as it would then be 
possible to see how the vortices respond to different types 
of pinning site. This type of experiment has been under- 
taken by Beleggia et al^ but so far the only coinparisons 
with theory have been visual. As shown in rcL [l^, much 
of the information on the structure of the vortex is con- 
tained in the contrast of the image rather than its visual 
appearance. Here we use images of flux vortices taken 
from MgB2 to make a quantitative comparison with sim- 
ulations to assess how much information can be obtained 
using transmission electron microscopy. We first compare 
simulations of defocussed images with those obtained ex- 
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FIG. 1. (a) Experimental arrangement for imaging flux vor- 
tices. The electrons are deflected by the component of the 
B-fleld from the vortices normal to the electron beam giving 
a black-white feature in an out-of-focus image, (b) The speci- 
men geometry. The MgB2 specimen was mounted to a copper 
post glued to a standard 3 mm diameter copper ring at an 
angle of 45°. An electron transparent window was then cut 
by focussed ion-beam milling. 



pcrimcntally and then discuss the use the transport of 
intensity equation which in principle allows the magnetic 
structure of flux vortices to be derived directly from the 
experimental images. 



II. MAGNESIUM DIBORIDE 

MgB2 was discovered to be a superconductor in 2001—. 
It has a transition temperature Tc = 39 K and a hexago- 
nal crystal structure (space group 191: PQ/mmm) com- 
posed of alternating layers of magnesium and boron with 
lattice parameters a ^ b ~ 3.086 A and c = 3.542 A. 
It is an unusual superconductor as it has two sources of 
electrons which contribute to superconductivity: one as- 
sociated with the a bonding from the boron p^y orbitals 
and the other associated with tt bonding from the boron 
Pz orbitals. As a result, the penetration depth varies 
with the applied field as this changes the proportion of 
£7 to TT electrons contributing to superconductivitjii^J^. 
In this experiment, the specimen was thinned to electron 
transparency in the c direction and since flux vortices 
penetrate parallel to the thin direction of the crystal, we 
make a low-field (~ 30 Oe) measurement of the penetra- 
tion depth in the a6-plane, Ao^. 

Only a few papers give values for the low-field penetra- 
tion depth for MgB2 resolved in crystallographic direc- 
tions. Cubitt et alr^ use small-angle neutron scattering 
in which measurements were extrapolated to zero field 
(the lowest was made at 0.1 T) to give Aab = 82 ± 2 nm 
at 2 K. Very low field (~ 1 ^T) radio-frequency mea- 
surements were also made by Manzano et ali^ giving 
Aab = 110-130 nm and Ac = 210-280 nm and a coher- 
ence length ^(j6 = 5.5 nm as the temperature approaches 



absolute zero. They also show that that Aab increases by 
8.0± 1.6 nm as the temperature is increased from 1.35 K 
to 10.8 K, the temperature at which the experiments were 
conducted here. 



III. METHODS FOR IMAGING FLUX 
VORTICES USING ELECTRON MICROSCOPY 

As flux vortices arc magnetic objects, they deflect but 
do not absorb the electron beam from the transmission 
electron microscope and so affect only the phase and not 
the intensity of the electron wavefunction. Consequently 
they are invisible in an in-focus image. There are two 
main imaging modes by which flux vortices can be visu- 
alised, giving access to the phase information. The first 
is off-axis holographj*i^ where a positively charged wire 
(called an electron biprism) is placed beneath the speci- 
men and used to interfere electrons which passed through 
the specimen with those which passed through vacuum to 
produce an interference pattern called a hologram. From 
this, the phase of the electron wavefunction can be cal- 
culated directly and differentiating the phase gives the 
B-field or more precisely, the induction-thickness prod- 
uct: the component of the B- field normal to the electron 
beam integrated along the path of the beam. 

The other method is out-of-focus (also called Fresnel) 
imaging where images are taken out of focus and the 
vortices appear as black-white features^. These can be 
compared with simulations or, if several images are taken 
at different defoci, the phase shift can be reconstructed 
using the transport of intensity equation but this is not 
as direct as holography as the intensity in a out-of-focus 
image is not sensitive to the phase or the phase gradi- 
ent but only to the curvature and higher derivatives of 
the phase. As the induction-thickness product is propor- 
tional to the first derivative of the phase, this has the 
consequence that it can only be determined to within an 
additive constant. It might be thought, therefore, that 
off-axis holography would be the better method for inves- 
tigating the structure of flux vortices as it gives a direct 
measurement of the phase shift but in fact, it has only 
been used to examine vortices in niobium by Bonevich et 

The main reason why out-of-focus imaging is the pre- 
ferred method for investigating flux vortices lies in the 
way Shot noise in the image is transferred to noise in the 
recovered phase. In the Supplemental Informatioi>ii we 
show that if the m-th derivative of the phase is denoted 
0'^'") then the noise associated with the phase derivatives 
recovered by holography is 

I /2m+l 

where N is the number of counts per reconstructed 
pixel, V is the visibility of the holographic fringes (the 
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difference between the maximum and minimum intensity 
divided by the sum) and X is the size of one pixel. 

On the other hand, the noise associated with the phase 
derivatives obtained from Fresnel images is 

Using reasonable values oi V = 30%, N = 200, A = 
0.00197 nm (the wavelength of 300 kV electrons) gives 
the noise in the curvature of the phase recovered by off- 
axis holography as cr^j^' ^ = 2.3 mrad/nm^ for the pixel 
size X = 16.9 nm used in this experiment. The noise can 
be reduced by increasing the pixel size and if X = 50 nm 
is used for the same number of electrons arriving at the 
detector, the noise is cr^j' ^ = 0.09 mrad/nm^. This is 
close to the lowest resolution at which individual vortices 
can be resolved as the penetration depth of niobium is 
52 nm. 

Simulations for a 250 nm thick specimen tilted to a = 
50° give a maximum phase-curvature of 0.47 mrad/nm^ 
for niobium (taking A = 52 nm [l^) and 0.12 mrad/nm^ 
for MgB2 (taking A = 110 nm) so it is clear that Bonevich 
et al. were working near to the limit of what is possible 
using electron holography, as they acknowledge^. We 
acquired electron holograms from MgB2 and found that 
although the average B-field from the array of flux vor- 
tices could be observed, individual vortices could not be 
identified. 

On the other hand, the noise expected in a phase 
reconstruction from defocussed images is a^/^ ^ = 

0.03 mrad/nm^ for A/ = 1.1 cm and X = 16.9 nm 
(used in this experiment) which is an improvement of 
2 orders of magnitude compared with holography at the 
same resolution. This shows why out-of-focus imaging is 
the more commonly used method for imaging flux vor- 
tices. This mode of imaging is often referred to as 'semi- 
quantitative' but, as we show here, provided that the im- 
ages are energy-filtered and recorded on a medium with a 
linear response, they contain all the information needed 
to make a quantitative comparison with simulations. 



'Lorentz' lens and a Gatan imaging filter. The micro- 
scope was operated in low-magnification mode at a nom- 
inal magnification of 105 x with the main objective lens 
turned off and the image was focussed with the diffrac- 
tion lens. This proved more convenient than imaging 
with the Lorentz lens as a wider range of defoci could 
be accessed and in this mode, the selected area aper- 
tures become objective apertures and vice versa whereas 
in Lorentz mode, neither set of apertures is in the cor- 
rect plane. The images were energy filtered so that only 
electrons which had lost 0-20 eV on passing through the 
specimen contributed to the image and an aperture was 
used so that only the 000 beam and the low-angle scat- 
tering from the vortices contributed to the image and 
the other crystallographic beams were excluded. The 
sample was cooled using a Gatan liquid-helium cooled 
'IKHCHDT3010-Speciar tilt-rotate holder which has a 
base temperature of 10 K. 

The defocus and magnification were calibrated by ac- 
quiring images with the same lens settings as the original 
images from Agar Scientific's 'S106' calibration specimen 
which consists of lines spaced by 463 nm ruled on an 
amorphous film. The defocus was found by taking dig- 
ital Fourier transforms of the images acquired from the 
calibration specimen and measuring the radii of the dark 
rings which result from the contrast transfer functioii^S. 

A thickness map of the specimen was created by divid- 
ing an unfiltered image by an energy-filtered image and 
taking the natural logarithm^! which gives the thickness 
parallel to the electron beam, as a multiple of the in- 
elastic mean free path, A^. To determine A^, an electron 
hologram was taken at room temperature at an edge of 
the specimen which gives a phase shift proportional to 
the thickness, = CeVqI. Ge is a constant which de- 
pends only on the microscope voltage and has the value 
6.523 X 10^ mr^Y-^ at 300 kV. Vb, the mean inner po- 
tential, was calculated as Vq = 17.71 V from theoretical 
scattering factors given in ref. [22l. giving A^ = 152 ±2 nm 
and the thickness, varied from 330-400 nm across the 
field of view of Fig. [2][a). (Ideally the thickness of the 
whole specimen would have been determined by electron 
holography but the field of view of the interference region 
was not sufficiently large.) 



IV. EXPERIMENTAL METHODS 

MgB2 single crystals were synthesised by Dr J. Karpin- 
ski via the peritectic decomposition of MgNBg as de- 
scribed in ref. [l9l The samples were thinned to 250 nm 
in the c-direction so that they were electron transpar- 
ent using a Helios Nanolab dual-beam focussed ion beam 
microscope. To tilt the sample to a high angle [a in 
Fig. [TJa)), the specimen was mounted on a tilted cop- 
per post attached to a standard 3 mm copper ring as 
illustrated in Fig. H^b). 

Microscopy was undertaken using a Philips CM300 
transmission electron microscope operated at 300 kV 
equipped with an electron biprism for holography, a 



V. COMPARISON OF DEFOCUSSED IMAGES 
WITH SIMULATIONS 

We first compare defocussed images with simulations 
based on the London model^. of flux vortices. The London 
model describes the B-field associated with a fiux vortex 
but gives an unphysical divergence of the B-field as the 
centre of the vortex is approached. This is corrected by 
the Clem model^^ where the B-field profile is rounded 
near the centre of the vortex in a region called the vortex 
core which has a radius given by the coherence length, 
^. The London model should be a good approximation 
for MgB2, however, as S^ab is about 6 nm whereas the 



penetration depth is about 100 nm and the pixel size 
of the images acquired here was 16.9 nm. 

The phase shift, (j){x,y), an electron experiences on 
passing through a London vortex has been calculated by 
Beleggia and Pozzi^i using a coordinate system in which 
the electron beam approaches the specimen in the z di- 
rection and x,y denote the image coordinates. It should 
be noted that their equation correctly accounts for the 
boundary conditions on the B-field at the interface be- 
tween the superconductor and the vacuum so that the 
spreading of the B-field near the interface and the effect 
of the stray field are both contained in the model. If the 
intensity in the in-focus image is uniform over the field 
of view (a good approximation here since the in-focus in- 
tensity varies by about 5% over the size of each vortex 
image whereas the out-of-focus intensity varies by about 
40% at the lowest defocus used here as can be seen from 
Fig.HJc)), the in-focus electron wavefunction can be writ- 
ten as ip = e**^ and from this, any image which could be 
taken using the electron microscope can be simulated. 

A defocussed image can be simulated by using the 
Fresnel-Kirchoff integral^! to propagate the in-focus elec- 
tron wavefunction by a distance A/, called the defocus. 
This is most easily accomplished by Fourier transform- 
ing the in-focus wavefunction in x and y but not z and 
multiplying by the phase factor e""^^/^'' (where A is 
the electron wavelength and k the spatial frequency) , in- 
verse transforming and taking the square modulus. (Here 
and throughout we use the Fourier transform convention 
f{k) = f{x)e-^-^''''=dx.) It should be noted that 
under the conditions which vortices are imaged, the con- 
tribution from the spherical aberration of the lens is neg- 
ligible. Here the diffraction lens was used to focus the 
image which has a spherical aberration coefficient of sev- 
eral metres but the coefficient would need to be ~ 10^ m 
before the aberration and defocussing terms were com- 
parable. 

The variable parameters used to model the images 
were: the penetration depth Aab, the position of the vor- 
tex in the image x, y, the thickness of superconducting 
material i, the tilting angle of the specimen a and the 
angle 9 of the axis of the vortex as seen in the image. 
These were chosen to minimise the sum-square difference 
between the experimental and simulated images using 
the simplex algorithm given in ref. [25l . A similar method 
was used to investigate the properties of pn junctions in 
semiconductors by Twitchett et al^. 

The position of each vortex x, y and the rotation an- 
gle of the vortices 6 given by the simplex algorithm are 
readily checked from the images. To give an independent 
measure of a, the angle between the plane of the copper 
ring and the specimen was measured in the focussed-ion 
beam microscope as 55 ± 1°. However, in the electron 
microscope, the specimen was in shadow when the ring 
was horizontal and a tilt of —8.85° was required to re- 
veal the electron transparent window. This would give 
a = 46° if the axes of rotation were the same but since 
it is likely they were different by several degrees, a will 
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FIG. 2. (a) Defocus series taken in a field of 27 G at 10.8 K. 
The defoci are (top to bottom) A/ = 2.24, 1.10, -1.10 and - 
2.07 cm. (b) The same defocus series with simulated vortices 
subtracted for the case where all the variable parameters are 
optimised. The other fits mentioned in the text give very 
similar images, (c) Linescans at each defocus taken perpen- 
dicular to the axis of the vortex indicted by the white box in 
(a). FT: fitted thickness - simulation in which the thickness 
was used as a variable parameter giving / = 543 nm, a = 58° , 
Aab ~ 126 nm. CT: calibrated thickness - simulation in which 
the thickness was fixed at its calibrated value, I = 372 nm, 
giving Q = 51°, Aab = 107 nm. DL: dead layer - simulation 
for which the thickness was fixed at I = 300 nm to mimic a 
dead-layer of non-superconducting material, giving a = 49°, 
Aab ~ 97 nm. E - experimental data. 



be larger and we estimate a = 50 ± 5° . 

The specimen thickness parallel to the electron beam, 
I, was measured to ±2 nm using a combination of elec- 
tron holography and energy-filtered imaging as explained 
in section HVl The thickness of the specimen parallel to 
the thin direction, t, is related to the thickness parallel to 
the electron beam, /, via t = I cos a and the thicknesses 
quoted here are those parallel to the electron beam unless 
otherwise stated. This method gives the total thickness 
of the specimen but does not account for the possibility of 
'dead-layers' of non-superconducting material produced 
by ion-beam thinning at the surface of the specimen. 
This is the largest of the uncertainties in this experiment 
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FIG. 3. Scatter graphs showing the results of the parame- 
ter optimisation generated by the simplex algorithm. In (a), 
the specimen thickness parallel to the electron beam, I, was 
treated as a variable parameter. In (b), it was fixed at its 
calibrated value. In (c), the thickness was reduced by 50 nm 
(72 nm parallel to the electron beam) from its calibrated value 
to simulate a 'dead-layer' of non-superconducting material. 
The horizontal line shows the mean and points within the 
shaded region are 1 standard deviation from the mean. 



and its effect is discussed later. 

We began by fitting all the parameters to the experi- 
mental data. Fig. ^a.) shows the defocus series and the 
black-white circular objects are the flux vortices. The 
diagonal stripes in the image arc thickness undulations 
unintentionally introduced by focussed ion beam milling. 
As discussed in ref. which used the same sample, they 
represent thickness undulations of about 1 nm. Vortices 
which were strongly affected by the contrast generated by 
these undulations were excluded from the fit. (b) shows 
the same images but where the fitted vortices have been 
subtracted. The fact that the vortices have been very ef- 
fectively removed indicates that the fit is good, (c) shows 
linescans at each defocus taken perpendicular to the axis 
of a representative vortex indicated in (a) by the white 
box and it can be seen that there is a good match between 
simulation and experiment. 

Fig. ^a) shows scatter graphs of the optimal values 
of the parameters for each vortex. The mean penetra- 
tion depth is Aab = 127 nm and the standard deviation 
is 12 nm. There are no obvious outliers in the scatter 
graph so assuming that all the vortices have the same 
penetration depth gives Aab ~ 127 ± 3 nm. 

To our surprise, the simplex algorithm gave thick- 
nesses of superconducting material which were consid- 



erably larger than the calibrated values. The average 
calibrated thickness was 365 nm but the average fitted 
thickness was 498 nm. We were expecting the fitted 
thicknesses to be shorter as there may have been damage 
introduced by the focussed ion beam milling resulting in 
non-superconducting material on the specimen surfaces. 
The problem seems to arise because the specimen tilt 
and the specimen thickness are to some extent comple- 
mentary: if a vortex has a given projected length normal 
to the electron beam, this could be the result either of a 
thick specimen with a low tilt angle, a, or a thin spec- 
imen with a high tilt angle. The algorithm would then 
adjust the penetration depth to account for this. 

To assess the extent to which a change in thickness 
changes the fitted penetration depth, we fixed the thick- 
nesses at their calibrated values and re-ran the simplex 
algorithm. This resulted in Fig. El^b) which gave a mean 
penetration depth of Aab = 113 nm and a standard devi- 
ation of 8 nm so that Aab = 113 ± 2 nm. Fig. ^c) shows 
that this change has almost no effect on the quality of 
the fit. 

Finally we assess the effect of a dead-layer of non- 
superconducting material on the surfaces of the sample. 
In addition to measuring the total sample thickness via 
a combination of energy-filtered imaging and holography 
as described in section HVl we also took convergent-beam 
diffraction patterns from 5 regions of the sample under 
two-beam conditions. These were used to find the thick- 
ness of the cr ysta lline component of the sample as de- 
scribed in ref. |20| . The difference between the two mea- 
surements gives the thickness of any amorphous layers 
present but we found that to within the experimental 
error of ±10 nm, there was no difference. 

We repeated the simplex procedure with the supercon- 
ducting thicknesses fixed at 50 nm (72 nm parallel to the 
electron beam) below their calibrated values which we 
consider to be the largest plausible dead-layer thickness. 
Again, Fig. [2jc) shows that the fit is still good and this 
yielded an average penetration depth of Aab ~ 101 nm 
with a standard deviation of 8 nm. Averaging over all 
the vortices gives Aab = 101 ± 2 nm. Thus we can say 
that at 10.8 K, Aab lies between 100-130 nm and is more 
likely to be between 100-115 nm. 



VI. DERIVATION OF THE VORTEX 
STRUCTURE USING THE TRANSPORT OF 
INTENSITY EQUATION 

In the previous section it was shown that simulated 
images based on the London model of flux vortices in 
MgB2 provided a good fit to the experimental images. It 
would be advantageous, however, to derive the phase and 
the B-field from the images directly so that no aspects 
of the model are assumed. In principle, this can be done 
by reconstructing the phase from the defocussed images 
using the transport of intensity equation^!. 

The transport of intensity equation is a re-expression 
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of the Schrodingcr equation in terms of the intensity, /, 
and phase, 0, of the electrons combined with the con- 
dition that there is a constant flow of electrons. Using 
the same coordinate system as section IVl where the elec- 
tron beam approaches the specimen in the z-direction, 
has wavelength A and the image (located at z = 0) has 
coordinates x, y, the transport of intensity equation is: 

V.,.(/V.,</.) = -(^^)|^ (3) 

When the in- focus image has a uniform intensity, Iq, 
(a good approximation for the images acquired here) the 
equation can be simplified to Poisson's equation: 



2tt 



dz 



(4) 



It can then be solved Fourier transforming in x and y 
but not z to give the Fourier transform of the phase as: 



1 1 
27rA/o fc2 



F.T. 



5/ 
dz 



(5) 



where F.T. stands for the Fourier transform of the term 
in square brackets. The rate of change of the intensity in 
the electron beam at the exit-plane of the specimen z = 
may be estimated by taking two images with equal and 
opposite defocus, ±A/, so that || ~ ^^^^ 2"a/~'^'^^ • "^^"^ 
phase in real-space can then be found by calculating the 
inverse transform of Eqn. [5] numerically. The singularity 
at fc = in Eqn. [5] is indicative of the fact that the 
zero of phase is arbitrary so a constant can always be 
added without any physical consequences. Here the pixel 
at fc = was set to zero so the average phase in the 
image was zero. It can also be seen from Eqn. 2] that 
the recovered phase is not just insensitive to an additive 
constant but also to the addition of a phase ramp. 

Ideally one would want to measure the B-field from 
the vortices. The quantity closest to this which can 
be derived from electron micrographs is the induction- 
thickness product, '(Bt)(a;, y)': the component of the B- 
field normal to the electron beam integrated along the 
length of the beam. It is related to the derivative of the 
phase via: 



(Bt)(a;,y) 



Bx 



dz 



h 
2tw 



~d(t)/dy \ 
d4>/dx J 



(6) 



For a conventional magnetic material where the B-field 
is constant through the thickness of the specimen and the 
stray field is negligible, the induction-thickness product 
is simply the component of the B-field normal to the elec- 
tron beam multiplied by the thickness of the specimen. 
In the case of flux vortices, the stray field is not negligible 
and is included in the induction-thickness product. As 
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FIG. 4. (a) Original image taken with A/ = 1.10±0.02cm.(b) 
Original image taken with A/ = —1.10 ± 0.02cm. (c) The 
magnitude of the induction-thickness product derived using 
the transport of intensity equation, (d) Simulated induction- 
thickness product for Aab = 113 nm using the calibrated 
thicknesses. 



mentioned above, the transport of intensity equation is 
insensitive to the addition of a phase ramp so the induc- 
tion thickness product can only be recovered to within an 
additive constant. As the simulations were carried out 
with the same boundary conditions as the reconstruc- 
tions, both have the same constant offset. 

Fig. m shows two images, equally disposed either side of 
focus with A/ — ±1.10 cm. A longer defocus series was 
obtained but the transport of intensity equation works 
best for images closest to focus due to the approximation 
to the gradient of the intensity with defocus explained 
earlier. From these images, the modulus of the induction- 
thickness product shown in (c) was derived, (d) shows 
a simulation of the induction-thickness product for the 
same vortex array using Aab = 113 nm and it can be 
seen that it is much more strongly peaked at the centre 
of the vortices. This could either indicate that the core 
of the vortices was much wider than previously thought 
or that the transport of intensity equation smooths the 
B-field profile. 

To ascertain which of these possibilities was correct, 
we used the simulated phase shift to produce simulated 
defocusscd images and then applied the transport of in- 
tensity equation to these. Fig. [SJa) shows simulated pro- 
files of the induction-thickness product across the cen- 
tre of a vortex and it can be seen that it is indeed the 
use of the transport of intensity equation which broadens 
the B-field profiles. This effect is caused by subtracting 
two images close to focus to approximate the gradient 

9| ~ ^^^'^''2Af^^^^ ^^'^ ^^"-"^ ^ ^^^^ '^^ coherence in the 
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FIG. 5. (a) The modulus of the induction-thickness product 
{Bt) for a simulated London vortex (solid line) and recon- 
structions based on pairs of simulated images at the defoci 
shown. A/ — ±1.10 cm was used in this experiment, (b) 
Comparison of the induction-thickness product (Bt) recon- 
structed from experimental data for a single vortex with sim- 
ulations for different values of Aab- (c) Induction-thickness 
profile derived from experimental data using the transport of 
intensity equation together with simulated profiles. In the 
first, the thickness parallel to the electron beam was fitted as 
part of the simplex routine yielding I = 530 nm, a = 52°, 
Aab = 123 nm. In the second simulation the thickness was 
fixed at the calibrated value of I — 373 nm, yielding a = 43°, 
Aab = 104 nm. 



electron beam as suggested in ref. |27|. 

Thus, the transport of intensity equation does not im- 
mediately yield the B-field profile of a flux vortex. To 
compare experiment and theory, it is first necessary to 
simulate defocussed images from the theory and then 
generate a magnetic profile by applying the transport of 
intensity equation to these. Fig. Eljb) shows that when 



this is done, a close match is obtained between the ex- 
perimental profile and simulations based on the London 
model and that if the penetration depth were the only 
quantity subject to uncertainty, it would be possible to 
determine a best fit to within ±5 nm. 

As before, however, there is an ambiguity between 
the specimen thickness and the specimen tilt angle and 
Fig.m^c) shows that very similar profiles can be obtained 
for different sets of parameters. The best values for the 
penetration depth were most likely found by direct com- 
parison between experimental and simulated images in 
section fVl rather than the more convoluted process of ap- 
plying the transport of intensity equation to experimental 
images and to simulated images and comparing the re- 
sults. If this is done, however, it yields A^^ = 125 ±4 nm 
when the thickness is treated as a variable parameter 
and Aab = 112 ± 3 nm when the thicknesses are fixed 
at their calibrated values. As expected, these values are 
very similar to those obtained by comparing the experi- 
mental images with simulations in section [V] 



VII. SUMMARY AND CONCLUSIONS 

Transmission electron microscopy gives a method to 
investigate the magnetic structure of individual fiux vor- 
tices in type-II superconductors at low magnetic field 
and allows absolute measurements of the penetration 
depth to be obtained from individual vortices irrespec- 
tive of whether they form a regular array. In section IIIII 
we showed that out-of-focus imaging is a more sensitive 
method for measuring the B-ficlds associated with fiux 
vortices compared with off-axis holography. 

We investigated the magnetic structure of flux vortices 
in MgB2 by comparing simulations based on the Lon- 
don model of fiux vortices with out-of-focus images and 
found these gave a good fit indicating that the London 
model gives a good description of fiux vortices in MgB2. 
The main uncertainty in the fits was the fraction of the 
sample thickness which was superconducting. Treating 
the thickness as a variable parameter in the fits and av- 
eraging over 17 vortices yielded a penetration depth of 
Aab = 127 ±3 nm although this gave an implausibly large 
specimen thickness: an average of 498 nm as opposed to 
the calibrated average of 365 nm. Fixing the thicknesses 
at their calibrated values gave Aab = 113 ± 2 nm and re- 
ducing the thicknesses from their calibrated thicknesses 
by 50 nm to mimic a dead layer of non-superconducting 
material gave Aab = 101 ± 2 nm with only minimal 
changes to the goodness of fit. Thus at 10.8 K and in 
a field of 27 Oe, Aab lies between 100-130 nm and is 
more likely to be between 100-115 nm. 

The most reliable meaurement of the penetration 
depth is likely to be the small-angle neutron scatter- 
ing experiments from single crystals by Cubitt et alt^ 
which gave Aab = 82 ± 2 nm at 2 K. According to the 
radio-frequency measurements of Manzano et alr^, this 
will increase to 90 ± 2 nm at 10.8 K, the temperature 
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at which the measurements were made here. Although 
radio-frequency measurements are usually used to derive 
changes in the penetration depth, Manzano et alr^ were 
able to estimate absolute values of Aab = 110-130 nm 
at 1.35 K which wiU increase to Aab = 118-138 nm at 
10.8 K using data from polycrystals and values for the 
anisotropy. Our values lie between these two measure- 
ments. 

We also investigated the use of the transport of in- 
tensity equation which in principle allows the magnetic 
fields in the specimen to be derived without reference 
to a model. We found that although this gave a useful 
overview of the B-fields from the vortices, the B-field pro- 
files are broadened due to the way the change of intensity 
with defocus is approximated. This must be taken into 
account when comparing the reconstructions with models 
so the technique does not provide a model-free method to 
determine the magnetic structure of vortices. It should 
be noted that if off-axis holography could be used, the 
magnetic profile would also be broadened due to the size 
of the Fourier mask used in the reconstruction process 
and this too would need to be accounted for. 

The coherence length S,ab in MgB2 is only 5.5 nm [l3| 
and we did not observe any rounding of the B-field profile 
which could be ascribed to the coherence length. For 
a material with a larger coherence length like niobium, 
where the coherence length is 39 nm |18l |. observing the 
effect on the B-field profile should be quite feasible. 

Quantitative measurements of the magnetic structure 
of flux vortices is only one aspect of the information 
which can be obtained using this technique. Informa- 
tion on the pinning landscape of the specimen can be 
derived from measurements of vortex motion^ and we 
are currently preparing a paper on this subject. 
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Fourier transforms referred to are taken in x and y but 
not z. 



A. Noise in Phase Reconstructions from Holograms 

Here we show how noise present in an off-axis hologram 
is transferred to the reconstructed phase. The standard 
expression's for the intensity in a hologram from a pure 



phase object is given by 



lix) = 1 + y cos {2ttQx + 0(a;)) + n{x) 



(7) 



where V is the visibility of the holographic fringes, Q 
is the carrier frequency, (f>{x) is the true phase shift (as 
opposed to that measured experimentally) and we have 
added noise, n{x) which is the noise-to-signal ratio for 
each pixel (note that the average intensity per pixel has 
been normalised to 1). We now follow the usual proce- 
dure to reconstruct the phase from the hologram. 

The first step is to take the Fourier transform of the 
intensity of the hologram which gives: 

I{k) ^ 5{k) + (^F.T.[e**(")] * S{k - Q) 

+ F.T.[e-**(")] *(5(fc + g)^ +n(fc) (8) 

where F.T. denotes the Fourier transform of the term 
in square brackets, * is a convolution and n{k) is the 
Fourier transform of the noise. It can be seen that the 
Fourier transform is composed of a central peak and two 
sidebands centred at k = zLQ. 

Next, the origin of the transform is shifted to the centre 
of the sideband at fc = Q by convolving both sides of the 
above equation by S(k+Q). The sideband is then isolated 
by multiplying by a top hat function h(k) which has the 
value 1 over the size of the sideband and is zero elsewhere 
to give a new function J(fc): 



J(fc) = ( ^F.T. 



n{k) * 5{k + Q) h{k) (9) 



The inverse transform gives: 



VIII. SUPPLEMENTAL INFORMATION - THE 
EFFECT OF NOISE ON THE RECONSTRUCTED 
PHASE 

Here we investigate how the noise in an image is trans- 
ferred to the phase recovered from these images. The 
coordinate system is oriented so that the electron beam 
travels along z and the image lies in the ly-plane. For 
simplicity, we consider 1-dimcnsional images although 
the results can be readily extended to 2-dimensions. The 



J{x) ^ (^-e^'^'(^) + n(a;)e2'^*'5^ j * h{x) (10) 

It can be seen that in the absence of noise, J(x) is 
proportional to the wavefunction of the electron beam as 
it exits the specimen, e**^. The inverse transform of the 
top- hat function, h[x), acts as a smoothing function and 
it controls the resolution of the reconstructed phase. Re- 
ferring to the Argand diagram in Fig. [BJ it can be seen 
that the phase measured by holography at a particular 
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B. Noise in Phase Reconstructions from Fresnel 
Images 

Here we show how noise present in pairs of defocussed 
images is transferred to the phase recovered using the 
transport of intensity equation. The transport of inten- 
sity of equation is: 



xy 



2tt\ dl 



(15) 



For a phase object with unit average intensity, this 
simphfies to 



FIG. 6. Argand diagram showing the how noise in a hologram, 
n{x), leads to noise in the recovered phase, 0. 



27r\ dl 



(16) 



position, X, is the sum of the true phase, 4> and an addi- 
tional phase due to the noise Acj) so that i^hoio = + A0. 
The geometry of Fig. |6] shows that provided the noise- 
to-signal ratio is much smaller than the fringe visibility, 
n/V 1, the additional phase shift is 



A(t>{x) =2'^sin{2TrQx-(f>{x)) (11) 

The noise in the phase recovered by holography is then 
the root-mean-square value of A0(a;) averaged over the 
image. Denoting the average over the image by (), the 
noise in the phase image is 



The derivative is found from two images equally dis- 
posed either side of focus so that 



f27r\ /(A/)-J(-A/) 
V.,</> = -(^-j (17) 

If / denotes the ideal intensity of the image and /cxpt is 
the intensity measured experimentally subject to a noise- 
to-signal ratio at each pixel of n{x, y), we have 



/expt(A/)=/(A/)+n(A/) 



(18) 



and 



y((A0? 



y2 



(n^x)) 



(12) 



The last expression applies if the noise is Shot noise 
and there are N counts per pixel. It is the same as the 
formula given by Lichte et al^ although obtained by a 
different method. 

Taking the m-th derivative (denoted by superscripts in 
brackets) of Eqn. [TT]in the approximation that there are 
many pixels for each holographic fringe, X <C 1/Q gives 



n^^Ux) 
2 ^^^sin(27rgx 



<P{x)) (13) 



The first derivative of the noise means the difference 
in the noise between two neighbouring pixels divided by 
the size of the pixel, X: n'{x) = {ti{x + X) — n{x)) /X. 
Subtracting the two noise terms from adjacent pixels in- 
creases the average noise by a/2 and so the noise in the 
?7i-th derivative of the phase due to Shot noise in the 
original hologram is: 



X' 



(14) 



/expt(-A/) = J(-A/) 



-A/) 



(19) 



Giving 



t^Frcsncl 



AA/ 



(n(A/) - n(-A/)) (20) 



As in the previous section, subtracting the two noise 
terms increases the average noise by a factor of -^2 so 
the noise in the m-th derivative of the phase recovered 
by the transport of intensity equation from images with 
an average count of N electrons per pixel and subject to 
Shot noise is 



X" 



'-XAf 



N 



(21) 
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